Burden of COVID-19 Risk in the Municipalities of East Java Province, Indonesia, and its Associated Factors: Spatial Analysis Approach

East Java Province has the fourth-highest number of COVID-19 cases among all other provinces Indonesia. This study aimed to examine spatial effect on confirmed and risk factors. Data were analyzed using Geoda software obtain Global Moran’s Index Local Spatial Autocorrelation (LISA) QGIS 2.8.1 make a map. I scatter plots also used exploring bivariate association between potential predictors. The statistics value shows clustering in across municipalities (Moran’s I=0.3986). A positive autocorrelation was observed population density = 0.2059), vaccination coverage 0.322), laboratories 0.2322), ratio health worker 0.1617), household 0.0866). In comparison, negative correlation Enforcement Restrictions Community Activities’ levels -0,2420), average family member 0.0115). LISA cluster map that there 3 hot spots (Surabaya, Gresik, Sidoarjo) cold (Sampang, Pamekasan, Sumenep).